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Abstract — This  paper  investigates  one  of  several  possible  means  of  deflection  of  a  large  near-Earth 
object  (NEO)  on  a  potential  collision  course  with  Earth  by  means  of  modest  velocity  changes 
(AK  »  1-200  cm/s)  applied  to  the  NEO  at  perihelion.  Given  the  present  uncertainty  in  the  geometry, 
topography  and  materials  properties  of  the  asteroids  and/or  comets,  a  deflection  mechanism  is  required 
which  is  relatively  insensitive  to  such  details.  Stand-off  nuclear  devices  present  an  option  to  fulfill 
this  criterion.  Momentum  can  be  transferred  by  several  mechanisms:  directly  by  means  of  the  kinetic 
energy  of  the  debris  or,  more  efficiently,  by  ablation  resulting  from  the  absorption  of  X-rays. 
Thermonuclear  devices  additionally  produce  high  energy  neutrons.  Using  approximate  models  of  the 
various  types  of  NEO  materials  and  tabulations  from  the  literature,  we  construct  effective  X-ray  and 
neutron  mass  absorption  coefficients.  The  latter  are  incorporated  into  a  simple  model  to  calculate  the 
impulse,  and  hence  the  velocity  change  imparted  by  the  various  nuclear  products  to  all  the  NEO  materials 
models. 


1.  INTRODUCTION 

Ever  since  1932  when  Reinmuth  discovered  1862 
Apollo,  the  first  Earth-crossing  asteroid  (ECA),  it  has 
been  understood  that  the  Earth’s  orbit  resides  within 
a  swarm  of  asteroids  and  comets.  It  was  also  realized 
that  the  orbits  of  these  Earth-crossing  asteroids  and 
comets  (ECCs)  can  be  perturbed  to  intersect  and  even 
impact  the  Earth’s  surface.  In  1971,  Gehrels  discov¬ 
ered  Daedalus  upon  initiating  a  systematic  wide  field 
search  for  ECAs.  Soon  after,  other  asteroid  surveys 
were  also  initiated  resulting  in  the  discovery  of 
numerous  ECAs  and  ECCs.  An  ECA  is  defined  [1,2] 
as  an  object  moving  on  a  trajectory  capable  of 
intersecting  the  Earth's  capture  cross-section  as  a 
result  of  on-going  (on  the  order  of  tens  of  thousands 
of  years)  long-range  gravitational  perturbations  due 
to  the  Earth  and  other  planets.  ECCs  may  be  defined 
in  an  analogous  way.  These  surveys  found  that  about 
90%  of  the  potential  Earth-crossing  (and  impacting) 
objects  are  considered  to  be  near-Earth  asteroids 
or  short-period  comets,  with  the  remaining  10% 
intermediate  or  long-period  comets  having  orbital 
periods  longer  than  20  years.  As  a  group  these  bodies 
are  called  NEOs  (near-Earth  objects),  which  can  be 
classified  as  having  orbits  that  closely  approach  and 
possibly  intersect  the  Earth’s  orbit  (~1  AU). 

A  model  of  the  true  population  of  ECAs  tested 
against  observed  discoveries  made  during  searches 
at  the  Palomar  46-cm  Schmidt  telescope  and  the 
University  of  Arizona  Spacewatch  telescope  [3] 
indicates  that  there  are  about  20  ECAs  larger  than 
5  km  in  diameter,  about  1 500  larger  than  1  km  in 
diameter  and  135,000  larger  than  100  m  in  diameter. 


There  are  also  thought  to  be  an  order  of  magnitude 
greater  number  of  smaller  NEOs  that  have  not  yet 
been  detected  [2,4].  While  none  of  the  approx.  1700 
suspected  ECAs  larger  than  1  km  in  diameter  cur¬ 
rently  presents  a  hazard,  they  are  nonetheless  capable 
of  developing  Earth-impacting  trajectories  over  a 
period  of  thousands  of  years.  This  population  model 
does  not  take  into  account  the  direct  hazards  from 
ECCs,  many  of  which  are  extremely  difficult  to 
detect.  These  numbers  present  but  part  of  the  poten¬ 
tial  hazard,  an  additional  source  of  Earth  orbit¬ 
crossing  objects  can  arise  from  debris  resulting  from 
tidal  break-up  (Roche  lobe)  of  both  long-  and  short- 
period  comets  in  Sun-grazing  orbits  [5].  This  latter 
class  of  NEOs  can  present  a  hazard  if  the  debris 
stream  intersects  the  Earth’s  orbit  at  one  of  its  nodes. 
The  potential  danger  to  the  Earth  from  comets  which 
undergo  chaotic  orbital  variations  represents  another 
facet  of  the  NEO  hazard  whose  statistics  have  yet  to 
be  understood. 

Despite  the  fact  that  NEOs  were  regarded  to  have 
collided  with  the  Earth  throughout  its  history,  up  to 
about  1980  there  had  been  little  concern  about  the 
likelihood  of  a  terrestrial  collision  with  an  ECA  or 
ECC.  However,  interest  in  ECAs  and  ECCs  was 
stimulated  largely  as  an  outgrowth  of  the  work  of 
Alvarez  el  al.  [6],  which  in  1980  proposed  that  an 
asteroid  or  comet  impacted  with  the  Earth  and 
initiated  a  global  disturbance  that  resulted  in  the 
mass  extinction  of  the  dinosaurs  as  well  as  other 
life-forms  over  a  broad  range  of  genera  and  species. 

The  recent  impact  of  comet  Shoemaker-Levy  9 
(SL-9)  on  Jupiter  has  underscored  the  fact  that 
near-Earth  objects  not  only  cross  planetary  orbits, 
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but  can  also  impact  them  with  devastating  conse¬ 
quences  [7].  The  SL-9  event  dramatically  presented 
the  possible  hazards  from  NEOs  in  a  context  hitherto 
regarded  as  unlikely.  To  this  end  it  is  prudent  to 
consider  analytical  models  which  address  orbit 
management  as  a  means  of  protecting  the  Earth 
against  a  direct  NEO  threat. 


2.  DEFLECTION  OF  LARGE  ASTEROIDS 

This  paper  addresses  one  method  of  deflecting 
large  asteroids  on  a  collision  course  with  the  Earth. 
For  example,  an  asteroid  of  around  20  m  in  diameter 
has  the  potential  destructive  capacity  of  over  1000 
MT  of  TNT.  Such  a  possibility,  though  remote,  is  not 
inconceivable,  and  has  been  the  topic  of  several 
recent  conferences  [8-1 3]  and  articles  [14-1 7],  Given 
sufficient  warning,  a  mission  could  be  launched  to 
apply  a  velocity  change  (A  V)  to  the  asteroid  in 
question  so  as  to  cause  it  to  miss  the  Earth.  Ideally, 
a  first-order  estimate  of  the  precise  value  of  the 
velocity  change  depends  on  the  asteroid’s  orbit, 
where  in  the  orbit  the  impulse  is  applied  and  the 
desired  miss  distance.  For  a  miss  distance  of 
10,000  km  and  interception  at  perihelion,  AFs  for 
typical,  known  asteroids  range  from  around  a 
centimeter  per  second  to  several  meters  per  sec¬ 
ond  [8,  18].  The  question  then,  is  how  best  to  apply 
the  required  impulse?  This  is  a  function  of  many 
parameters  which  include  the  mission  economics  and 
logistics,  the  technical  means  deemed  most  appropri¬ 
ate  and  the  materials  properties  of  the  NEO.  Subsur¬ 
face  momentum  coupling  from  conventional  (HE) 
and  nuclear  (NE)  explosives  has  been  discussed  by 
Remo  and  Sforza  [19],  while  a  more  general  review  of 
the  momentum  coupling  options  (pulsed  laser,  kinetic 
energy  impactors,  NE  and  HE  explosives)  can  be 
found  in  the  paper  of  Shafer  et  al.  [20].  The  mechan¬ 
ism  that  appears  to  be  the  least  sensitive  to  current 
uncertainties  in  composition  and  shape,  i.e.  to  the 
geology  and  materials  properties  of  asteroids,  is  a 
stand-off  nuclear  burst  [21,22].  Since  the  latter  irradi¬ 
ates  approximately  half  the  surface  area,  it  is  also 
less  liable  initially  to  add  additional  angular  momen¬ 
tum  to  the  asteroid  as  would  be  the  case  for  a 
misplaced  surface  burst,  penetration  device  or 
thruster.  It  can  be  shown  that  for  the  ratio  of  the 
rotational,  T,  to  the  kinetic  energy,  E,  to  be  small. 


that  TjE  =  (Fa/AF)2«  1;  FA  =  <u[/l]  where  o>  is  the 
angular  velocity  and  A  a  characteristic  length. 

A  simple  model  [22-25]  adapted  to  treat  stand-off 
momentum  transfer  is  discussed  in  more  detail  than 
previously  [22]  and  is  then  applied  to  specific  classes 
of  asteroid  analog  materials,  meteorites.  This  model 
is  not  meant  to  take  the  place  of  more  detailed 
calculations.  However,  the  latter  require  multi-group 
neutron  and  radiation  transport  and  equation  of  state 
data  coupled  to  a  full  hydrodynamic  code  in  order  to 
estimate  the  energy  deposition  and  impulse  given  to 
the  asteroid.  Since  accurate  models  for  asteroids  and 
materials’  properties  are  lacking  for  the  latter,  it  is  felt 
that  a  first-order  approach  is  preferable  at  this  time. 

The  energy  liberated  in  a  nuclear  explosion  in  space 
can  be  divided  into  that  coming  directly  from  the 
nuclear  reactions  and  that  from  the  heated  case  and 
weapon  material.  Table  1,  from  the  ICARUS  report 
[26]  and  based  on  the  information  in  Glasstone  [27], 
gives  the  fractional  distributions  of  the  energy.  Based 
on  the  nuclear  reaction  kinetics,  Glasstone’s  table 
applies  to  fission  devices,  although  this  is  not  explic¬ 
itly  stated.  The  prompt  gamma-rays  are  emitted  in 
the  first  instants  of  the  explosion,  the  delayed  gamma 
and  beta  radiation  come  from  the  fission  products. 
The  bulk  of  the  radiation  is  emitted  in  thermal 
X-rays.  For  a  thermonuclear  device  the  energy  release 
derives  approximately  equally  from  fission  and 
fusion  [27];  most  of  the  fusion  energy  is  in  the  form 
of  high-energy  neutrons.  The  neutrons  are  slowed 
down  by  collisions  with  the  fuel  and  case  and  can 
cause  fast  fission  reactions.  This  gives  rise  to  a 
spectrum  of  energy  values  of  the  surviving  neutrons 
which  is  not  given.  Below  we  estimate  the  impulse 
given  to  the  asteroid  by  the  various  sources  listed  in 
Table  1:  the  debris,  and  ablation  due  to  the  absorp¬ 
tion  of  gamma-rays.  X-rays  and  neutrons.  The  results 
are  also  presented  in  graphical  form  in  the  figures  for 
representative  values  of  the  opacities. 


3.  IMPULSE  DUE  TO  DEBRIS 

The  specific  impulse  from  the  debris  is  approxi¬ 
mately: 


where  MD  is  the  debris  mass,  L  the  distance  from  the 


Table  1(a).  Fractional  energy  distribution  (fission  device)  [26,27] 
Nuclear  radiation  Thermal  radiation 


Prompt  gamma-ray  0.001  X-rays  0.70 

Neutrons  0.01  u.v.,  visible,  i.r.  0.05 

Delayed  gamma-rays  0.02  Debris,  kinetic  energy  0.20 

Delayed  beta-rays _ 002 _ 


_ Table  1(b).  Fractional  energy  distribution  used  for  the  fusion  device 

Source _ X-rays _ Neutrons _ Gamma-rays _ Debris _ Other 

Fraction _ CL55 _ 020 _ O01 _ 020 _ 0.04 
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burst  to  the  asteroid  and  <Fd>1/2  is  the  RMS  velocity 
of  the  debris.  Here  and  in  what  follows  we  assume 
that  the  distance  from  the  burst  is  much  larger  than 
the  size  of  the  object.  The  RMS  velocity  is  related  to 
the  energy  in  the  debris,  ED,  by: 

\MDiVl')=ED  =  t]Y  (2) 

where  Y  is  the  total  yield  and  t?  the  fraction  of  the 
yield  appearing  in  the  debris.  The  specific  impulse 
also  is  given  by: 

ID  =  (Ma/A)AVd  (3) 


where  MA  is  the  asteroid’s  mass  and  A  its  (total)  area. 
The  velocity  change  is,  on  combining  eqns  (l)-(3), 
and  taking  r\  =  0.2: 


In  the  last  term  of  the  right  hand  side  of  the  above 
expression  the  yield  is  in  MT(TNT  equivalent)  and 
the  mass  is  in  kg. 

The  geometry  of  the  encounter  is  shown  in  Fig.  1 . 
The  fractional  solid  angle.  Aft  =  (A  /4FIL2),  depends 
on  the  ratio,  y  =  H/R,  of  the  stand-off  distance,  H, 
to  the  radius,  R,  for  a  spherical  asteroid,  and  L 
corresponds  to  the  tangential  distance  in  Fig.  1. 

With  reference  to  Fig.  1,  the  fraction  of  the  area 
irradiated  is  /=  1/2[1  —  cos  </>  —  y/(l  +  y)],  and 
Aft  =  //[y2  -(-  2y].  Similarly,  the  fraction  of  the  energy 
liberated  and  subtended  by  the  NEO  is 
g  =  1/2[1  -  sin  <t>  =  (y2  +  2y)1/2/(l  +  y)].  There  is 
thus  a  trade-off  between  a  more  uniform  irradiation 
and  the  energy  requirements.  For  example,  for  y  =  2, 
/  =  1/3  and  Aft  =1/24.  The  value,  y=v/2-l, 
(y  =  0.414)  for  which  /  =  0.15  =  Aft  and,  in  Fig.  1, 
(p  =  a  =  45°,  results  from  requiring  /  +  g  be  a  maxi¬ 
mum  (see  also  [18]).  Finally,  Hyde  [21],  using  a  simple 
model,  found,  upon  numerical  integration,  that 
y  =  0.5  close  to  the  above  geometric  optimum. 
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Taking  MD  =  300  kg  [23], 

A  VD  =  ZlLiliL  AO  /  T(MT)  cm/s.  (5) 

Ma 

Since  the  mass  of  a  kilometer  size  asteroid  can  exceed 
1012kg,  the  debris  contribution  can  be  neglected  on 
this  scale.  The  debris  can  be  effective  for  objects  of 
the  order  of  200  m  diameter  for  the  parameters  used 
in  this  example.  The  effectiveness  of  the  debris  in 
deflecting  large  diameter  objects  could  be  enhanced  if 
the  mass  of  the  debris  and  hence  the  coupling  to  it 
were  increased.  For  yields  in  excess  of  1  MT,  the  mass 
of  the  debris  presumably  increases,  assuming  it  is 
scaled  according  to  MD  ~  Y2  '  the  velocity  increment 
is  enhanced  over  that  in  eqn  (5)  only  by  AVDcc  K1 3. 


4.  IMPULSE  ANALYSIS  FOR  RADIATION 

We  first  describe  the  analytical  model  used  to 
estimate  the  impulse  [22], 

The  specific  impulse  (momentum  per  unit  area), 
/,  imparted  by  ablation  due  to  X-ray  or  neutron 
absorption  is  found  by  integrating  the  differential 
relationship: 

dl  =  pV  d.v  (6) 

where 

V2/2  =  t  (x)  U  (7) 

p  is  the  density,  e(.v)  and  ev  are  the  incident  energy  per 
unit  mass  and  vaporization  energy  per  unit  mass, 
respectively,  and  d.v  is  the  radial  differential  to  the 
NEO  surface.  e(x)  is  related  to  the  absorbed  fluence 
F,  and  the  opacity  p  (cm2/g),  by 

£=pF.  (8) 

Therefore 

7  =  \/  2  j"  '  (pF  —  tv)l  2p  d.v  (9) 

where  Xv  is  the  evaporation  depth. 

Additionally,  it  is  assumed  that  the  fluence  varies 
with  the  penetration  depth  according  to 

F  =  F0  exp(  -  ppx  /cos  9 )  (10) 

where  9  is  the  angle  between  the  incident  radiation 
and  the  normal  to  the  surface.  Assuming  a  nuclear 
device  explodes  at  a  distance  L  from  the  NEO’s 
surface,  the  incident  fluence  is, 

F0  =  rj  Y!4nL 2 

=  1.33  x  1 0S>7  Y (MT)AQ/£>2  J/cm2  (11) 
where  the  yield,  T,  is  expressed  in  MT  (TNT  equival¬ 
ent)  and  D,  the  asteroid  diameter,  in  km.  rj  denotes 
the  fraction  of  the  yield  in  X-rays,  neutrons  or 
gamma-rays.  The  temperature  variation  of  the 
opacity  is  ignored.  For  simplicity  the  obliquity 
factor,  cos  9,  is  set  equal  to  unity  which  is  valid  for 
not  too  large  values  of  the  angle  a  of  Fig.  1 .  For 


example,  using  x  =  0.414  and  a  =  22.5°,  instead  of  its 
maximum  value,  1  /cos  9  =  1.01. 

The  specific  impulse  can  be  easily  obtained  by 
using  the  scaled,  non-dimensional,  variables 

F*  =  F0/Fmm  ^  1 

Fmin  =  £v/p.  (12) 

Fmin  represents  the  minimum  fluence  that  is  required 
to  evaporate  the  material. 

The  evaporation  depth  is  given  by: 

Xy  =  (cos  9jpp)\n  FJ .  (13) 

With  this  form,  the  integral  can  be  readily  evaluated 
/*  =  2v/2{v/f'0*  -  1  -  tan-’I^/F?  -  U}-  (14) 

Note  that  the  material  properties  in  this  model  are 
entirely  contained  in  the  simple  ratios  of  eqn  (12). 

If,  instead  of  an  exponentially  decreasing  absorp¬ 
tion,  one  used  a  constant  absorption  model  up  to  X, 
and  zero  thereafter,  the  above  calculation  gives 

I*  =  J2 (F0*  -  1).  (15) 

These  and  other  similar  models  also  have  been  dis¬ 
cussed  by  Lawrence  [24,25].  In  the  asymptotic  regime, 
FJ»1,  and,  using  eqn  (12) 

/=2v/2F0/p  exponential  absorption  (16a) 
=  ^2F9/p  constant  absorption  (16b) 

t„,  does  not  appear  in  the  expressions  in  this  limit, 
but  of  course  must  be  taken  into  account  in 
determining  if  indeed  FJ  » 1 .  /  is  usually  measured 
in  taps  (g/cm-s). 

It  is  seen  that  in  the  asymptotic  limit,  other  than 
the  fact  the  specific  impulse  for  exponential  absorp¬ 
tion  is  twice  that  for  constant  absorption,  there  is  the 
same  parametric  dependence  on  the  fluence  and  the 
opacity  in  eqns  (16a)  and  (16b).  The  surface  materials 
properties  of  the  asteroid  are  expressed  through  p. 
The  fundamental  assumptions  of  the  model  are 
mono-dimensionality,  absorption  coefficients  inde¬ 
pendent  of  depth  (i.e.  uniformity  of  composition)  and 
temperature,  and  neglect  of  viscosity  and  thermal 
conduction. 

Hence,  for  a  given  yield,  energies  of  vaporization 
and  the  mass  absorption  coefficients,  the  impulse  can 
be  calculated  from  eqns  (12)  and  (14).  The  resulting 
velocity  change  of  the  asteroid  can  be  obtained  from 
the  relation, 

/  =  wAF  (17) 

where  m  is  the  areal  mass  density  of  the  asteroid.  In 
terms  of  Opik’s  [28]  shape  parameter,  5,  and  the  mass 
density,  p , 

m  =  pRjB.  (18) 

(5  =  3  for  a  sphere.) 
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Finally,  we  have 

AV  =  (2s/2/mtt)J7wyF*-\ 

-tan 'l[y/Fg  -  1]}.  (19) 
In  terms  of  the  non-dimensional  quantities, 

A  V*  s  A and  m*  =  mp,  (20) 

we  have, 

AK*  =  (2v/2/m*){v/F0*-l 

-tan'  Vf?  -1]}.  (21) 


As  discussed  in  the  APS  report  [23]  and  by 
Lawrence  [24],  a  useful  way  of  contrasting  experimen¬ 
tal  data  with  theory  is  to  plot  the  impulse  coupling 
coefficient, 


Cm  =  I/F0  (22) 

versus  the  incident  fluence,  F0.  The  units  of  Cm  are 
conventionally  taken  to  be  dyne-s/J.  The  impulse 
coupling  coefficient  has  long  been  used  extensively  to 
analyze  laser-induced  impulse  [24,29]  but  is  not  re¬ 
stricted  to  that  application.  In  terms  of  the  non-di¬ 
mensional  quantities  defined  above, 

CZ  =  I*/F$ 

AV*  =  (F*  /m*)CZ 

yfeCm.  (23) 

When  I*  is  given  by  eqn  (14),  C*  is  maximal  for 
F*  =  6.434,  when  C*  =  0.5124 [24],  For  the  problem 
under  discussion,  as  will  be  seen  below,  for  many  of 
the  values  of  the  parameters,  C*  is  not  optimal  in  this 
sense. 

In  the  next  sections  we  discuss  the  contributions 
to  the  opacities  due  to  gamma-rays.  X-rays  and 
neutrons.  Although  they  are  not  an  important 
contributor  for  our  purposes,  we  include  the  gamma- 
ray  contribution  for  completeness. 


5.  OPACITIES 


5.1.  Gamma -rav 

According  to  Glasstone  [27],  the  gamma-ray  spec¬ 
trum  from  a  burst  in  outer  space  arises  primarily 
from  the  delayed  gammas  produced  by  the  fission 
products,  the  prompt  gammas  being  absorbed  in 
the  case  and  nuclear  material  before  disassembly.  The 
delayed  gammas  have  an  energy  of  ~2MeV.  The 
opacity  used  in  eqn  (8)  is  defined  by:  p  =k/p,  where 
k  is  the  linear  absorption  coefficient  and  p  the  mass 
density.  Using  Glasstone’s  Table  8.75  [27]  and  density 
values  of  2.3,  7.9  and  11.3g/cm3  for  concrete,  iron 


Table  2.  Gamma-ray  mass  absorption  coefficient  (cmz/g)  [27] 
Energy  (MeV)  H2Q  Concrete  Iron  Lead 

2  0.049  0.050  0.042  0.046 

5  0.030  0.031  0.032  0.044 

10 _ 0.022  0.026  0.029  0.054 


Table  3.  X-ray  mass  absorption  coefficient 

_ (cmz/g)  [30,31] _ 

X-ray  energy  (keV) 


Al  191  25.5 

O  49  5.8 

Mg  158  20.4 

Si  244  33 

Fe  138  176 

Ca  637  96.5 

S  353  49.6 

Ni _ 180  215 


6.2 

3.2 

13.8 


and  lead,  respectively,  one  can  construct  a  table 
(Table  2)  for  py . 

The  values  in  Table  2  should  be  contrasted  with 
Glasstone’s  [27]  choice  (Section  8.79)  of  0.023  cm2/g 
as  a  value  to  use  for  most  materials.  For  the  estimates 
below  we  will  use  the  “universal”  value  of  0.05  cm2/g 
as  appropriate  for  a  2  MeV  gamma-ray. 

5.2.  X-rays 

We  will  assume  that  the  X-rays  are  emitted  at  a 
black-body  temperature,  T.  Glasstone  [27]  (Section 
7.8)  states  that  most  of  the  X-rays  are  emitted  in  the 
range  of  0.1-100  A  (0.12-120  keV),  the  temperature 
being  “several  tens  of  millions  of  degrees  Kelvin”. 
Expressing  T  in  units  of  keV,  the  black-body  maxi- 

Anax  =  (2.5/TkeV)  A 


Em  ax  =  5TleV  keV.  (24) 

The  ICARUS  report  [26]  used  £raax=  120  keV 
corresponding  to  T  =  24  keV,  but  this  appears  too 
high.  A  temperature,  T  =  2  keV  (2.3xl07K)  for 
which  £max  =  10  keV,  seems  more  reasonable,  at  least 
for  a  fission  device.  In  order  to  examine  the  sensitivity 
to  the  assumptions,  Table  3  lists  mass  absorption 
data  [30,31]  for  photon  energies  of  5,  10  and  25  keV; 
the  latter  value  scaled  from  that  at  10  keV  according 
to  px  ~  (l/£max)3.  The  choice  of  elements  is  discussed 
in  Section  6  below.  In  principle,  eqn  (9)  should  be 
averaged  over  the  Planck  distribution;  however,  these 
data  exist  in  tabular  rather  than  analytical  form  and 
it  was  felt  that  this  would  have  rendered  the  results 
less  transparent. 

Opacities  for  composite  materials  are  estimated  by 
summing  over  those  of  its  constituents  weighted 
according  to  their  atomic  percent,  a,: 

(25) 

The  X-ray  absorption  cross-sections  ( Q  =  pxA/N0, 
where  A  is  the  atomic  weight  and  N0  is  the 
Avogadro’s  number)  are  three  to  four  orders  of 
magnitude  greater  than  those  for  the  neutrons,  as  will 
be  discussed  below,  hence  their  penetration  depth  is 
correspondingly  less. 
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Table  4.  Neutron  total  cross-section 

_ (barns)  [33] _ 

Neutron  energy  (MeV) 
Element  5  14 

A!  2.4  1.7 

O  1.0  1.5 

Mg  2.0  1 .65 

Si  2.5  1.85 

Fe  3.7  2.5 

Ca  3.3  2.1 

S  2.6  1.8 

Ni _ 3/7 _ 2.7 


5.3.  Neutrons 

The  principal  reactions  involved  in  the  fusion  of 
the  hydrogen  isotopes  are  [32]: 

D  +  D-T(1.01  MeV)  +p  (3.02  MeV) 
->He3(0.82  MeV)  +  n(2.45  MeV) 

D  +  T-He4  +  n(14.1  MeV) 

D  +  He3-He4  +  p(14.7  MeV) 

T  +  T->He4  +  2n  +  1 1.3  MeV.  (26) 

The  two  D-D  reactions  occur  with  approximately 
equal  likelihood.  The  D-T  reaction  rate  for  tempera¬ 
tures  between  10  and  30keV  is  two  orders  of 
magnitude  greater  than  that  of  the  D-D  and  the  other 
reactions  and  thus  dominates  in  the  range  of  interest. 
In  what  follows  we  will  also  ignore  the  reactions 
between  the  by-products.  The  neutrons,  initially  at 
14  MeV,  undergo  collisions  in  the  D-T  fuel  (assumed 
to  be  half  consumed),  in  the  case  and  other  materials 
thus  acquire  an  energy  distribution  or  spectrum. 
Fortunately,  the  neutron  cross-section  is  relatively 
insensitive  both  to  the  energy  distribution  and  the 
material  in  our  range  of  interest.  For  the  fission 
neutron  spectrum  for  slowing  down  in  air  at  stp, 
Glasstone  [27]  (Fig.  8.96  and  Section  8.104)  uses 
a  constant  cross-section  of  2b  (\b  =  10~24 cm2). 
Assuming  this  applies  to  solid  density  matter 
(N  =  5  x  1022/cm3),  this  gives  a  mean-free-path  of 
10  cm.  For  the  more  energetic  neutrons  coming  from 
the  D-T  reaction,  taking  a  range  of  5-14  MeV  to 
partially  take  the  spectrum  into  account,  the  cross- 
sections  vary  up  and  down  by  no  more  than  a  factor 
of  two  from  the  above,  as  can  be  seen  in  Table  4  [33]. 
Again  the  choice  of  elements  is  dictated  by  the 


materials  of  interest.  For  the  estimates  below  we 
assume  therefore,  that  all  the  neutrons  are  14  MeV. 
The  D-D  neutrons  contribute  much  less  to  the 
overall  fusion  energy  due  to  their  lower  energy  and  to 
their  reduced  reaction  rate  and  probability  to  occur. 

Heating  by  neutrons  requires  a  more  detailed 
treatment  than  can  be  included  in  simple  models.  In 
particular,  the  elastic  scattering  of  the  neutrons  by  the 
ambient  ions  leads  to  recoil  of  the  latter.  These 
“knock-on”  ions  must  be  followed  in-turn  as  they 
slow  down  and  deposit  their  energy  in  the  bulk  of  the 
material.  The  net  result  for  the  specific  impulse  may 
be  intermediary  between  a  constant  absorption 
model  and  an  exponential  model.  Therefore,  both  the 
constant  absorption  and  exponential  model  results 
for  the  neutrons  are  contrasted  in  Section  7  below. 

6.  ATOMIC  COMPOSITIONS  OF  ASTEROID  MATERIALS’ 
ANALOGS 

The  asteroid  compositions  described  below  and  in 
Table  5,  will  be  used  to  compare  the  effects  of  the 
various  radiation  sources.  The  models  are  based  on 
meteorite  analogs  and  are  meant  to  be  illustrative 
and  are  tentative  pending  further  data  [34,35].  The 
chondrite  groups  [36],  are  subdivided  into  non- 
volatile-free  (NVF)  and  volatile  (S.  O,  C.  H)-free 
(VF)  in  order  to  account  for  the  possibility  of  out- 
gassing  of  a  surface  layer  during  the  first  of  several 
bursts. 

(1)  Iron-nickel 

Composition:  90%  Fe,  7-8%  Ni  (atomic  %) 
Density:  7.9g/cm3 

(2)  Enstatite 

Composition:  see  Table  5 
Density:  4g/cm3 

(3)  Ordinary  chondrite 

Composition:  see  Table  5 
Density:  3.4g/cm3 

(4)  Carbonaceous  chondrites 

Composition:  see  Table  5 
Density:  2. 1  g/cm3 

Opik’s(28]  values  for  the  vaporization  energies  of 
iron  and  stony  meteorites  are  8.26  x  1010  and 
8.08  x  10‘°ergs/g,  respectively,  while  Bronshten  [37] 
uses  the  value  of  £v  =  8.0  x  10‘°ergs/g  for  all  meteor¬ 
ites.  The  latter  value  will  be  adopted  here  as  well. 


Table  5.  Composition  (atomic  %)  of  chondrite  groups  [36] 


Element _ NVF _ VF 


A1  0.81  1.7 

Fe  16.44  34.6 

Mg  12.13  25.6 

Ca  0.81  1.7 

Ni  0.81  1.7 

S  5.12  — 

_ O _ 47.44 _ — 

NVF  =  non-volatile  free;  VF  =  volatile  free. 


NVF 


16.0 

1.2 

9.4 

14.9 

0.8 

1.7 

55.88 


VF _ NVF 


2.8  1.2 

22.0  11.39 

35.0  15.1 

1.7  0.99 

1.2  0.74 

1.48 

— _ 55.2 
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Table  6.  lOkeV  X-ray  and  14MeV  neutron  mass  absorption  co- 
_ efficients  (cm2/g) _ 


OO _ <F  N> 


Enstatite  NVF  44.9  4.8  x  10  2 

Enstatite  VF  83.5  3.6  x  10  2 

Ordinary  NVF  30.6  5.1  x  10  "2 

Ordinary  VF  63.4  3.8  x  10~2 

Carbonaceous  NVF  34.6  5.0x  10  2 

Carbonaceous  VF  70.8  3.8  x  10~2 

Fe/Ni  (90%/10%) _ 180 _ 2.8  x  IQ-2 


NVF  =  non-volatile  free;  VF  =  volatile  free. 


7.  RESULTS 


The  above  information  can  now  be  combined  into 
estimates  of  the  yields  required  for  given  velocity 
increments,  densities  and  diameter.  First  consider 
the  gamma-ray  contribution.  The  non-dimensional 
parameter  Fg ,  in  units  already  defined,  with  the  value 
of  tv  given  above,  and  on  using  eqns  (11)  and  (12),  is 
given  by: 


Using  the  parameters  for  gamma-rays,  F*  ~ 
(0.008)T(MT)/D:.  Thus  gamma-rays,  especially  for 
large  NEOs.  will  not  make  a  significant  contribution 
until  around  the  1000  MT  range. 

Equation  (21)  can  be  simplified  by  defining, 

y  =  m*AV*/2y/2 
<t>  =  (F*  -  1)1/2 

so  that 


(28) 


Here  y  is,  apart  from  the  numerical  factors,  the 
non-dimensional  momentum  transferred  to  the  NEO, 
while  <j>2+  1  is  the  non-dimensional  fluence. 
Expressed  in  terms  of  the  physical  quantities, 

y  =  ppDAV/l^/lt,  =  ppD  AV/4S.  (29) 


Fig.  3.  Specific  impulse  (ktaps)  vs  the  yield  parameter 
foAnr(MT)]  produced  by  lOkeV  X-ray  and  14MeV  ab¬ 
sorption  for  D  =  2  km  and  D  =  10  km  chondrite  and  Fe/Ni 
asteroid  models. 


Again  we  have  used  the  value  of  £v  given  at  the  end 
of  Section  6. 

The  X-ray  and  neutron  opacities,  averaged  over  the 
atomic  percents,  are  listed  in  Table  6.  Representative 
values  for  the  X-ray  opacities  are  35  cm2/g  for  non¬ 
volatile  free  chrondites,  about  twice  that  value  for  the 
volatile-free  material,  and  180  cm2/g  for  Fe/Ni  NEOs. 
The  corresponding  neutron  opacities  are  0.05  and 
0.028  cm2/g  for  all  chrondites  and  iron/nickel  NEOs, 
respectively.  From  eqn  (29)  the  density  and  opacity 
enter  as  a  product,  so  that  the  value  of  y  for  neutron 
interactions  is  approximately  the  same  for  the  chon¬ 
drites  and  Fe/Ni  NEOs. 

Since  y  is  proportional  to  the  opacity,  the  low 
opacities  for  neutrons  result  in  their  values  for  y 
falling  in  the  region  where  the  full  expression 
[eqn  (28)]  must  be  used.  The  X-ray  opacities  are 
sufficiently  high  that  their  values  for  the  non- 


Fig.  2.  Velocity  imparted  by  lOkeV  X-ray  and  14MeV  neutron  absorption  vs  the  yield  parameter 
[»/Afiy(MT)]  for  D  =2  km  and  D  =  10  km  chondrite  and  Fe/Ni  asteroid  models. 
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dimensional  momentum  parameter  fall  in  the  asymp¬ 
totic  region.  It  is  useful  to  present  this  information  in 
its  dimensional  form.  The  desired  velocity  increment, 
AV  (cm/s),  is  plotted  in  Fig.  2  versus  the  quantity 
tjy(MJ)AQ,  known  as  the  yield  parameter,  for  an 
X-ray  and  neutron  radiation  of  2  and  10  km  diameter 
chondrite  and  Fe/Ni  asteroids.  Figure  3  shows 
results  for  the  specific  impulse  for  the  same  par¬ 
ameters. 

8.  DISCUSSION  AND  CONCLUSION 

The  trends  seen  in  Figs  2  and  3  are  due  to  the  large 
differences  in  the  opacities  for  X-ray  and  neutron 
absorption.  This  behavior  can  be  explained  by  con¬ 
sidering  eqn  (28),  for  X-rays  y  » 1  and  <j)  «  y,  whereas 
for  neutron  absorption  y«  1  and  <p  ~  (3y)1,3.  In  this 
application.  X-ray  absorption  in  the  exponential 
model  can  be  approximated  by  its  asymptotic  limit 
whereas  neutron  effects  cannot.  In  the  asymptotic 
limit  both  the  velocity  change  and  the  specific  impulse 
are  inversely  proportional  to  the  inverse  square  root 
of  the  opacity,  proportional  to  the  square  root  of  the 
flux,  and  independent  of  the  heat  of  vaporization. 
The  effect  of  neutron  irradiation  requires  use  of  the 
full  expression  and  therefore  exhibits  the  steep  cut-off 
behavior  seen  in  the  figures,  which  occurs  when  the 
non-dimensional  flux  parameter  FJ  «  1 .  The  expo¬ 
nential  absorption  model  thus  gives,  for  X-ray  and 
neutron  absorption  respectively: 

F(MT)An/£>2  w  0.06(p»AF/48)Vx//C1 
«  1.2[1  +  (3y)2/J]/ijN 
«0.06[l+4y2]/i/N.  (30) 

The  last  line  of  eqn  (30)  is  the  result  for  the  constant 
absorption  model  for  neutron  absorption. 

As  discussed  in  Section  1,  the  contribution  to  AF 
from  the  debris  for  asteroid  diameters  greater  than 
1  km  can  be  neglected  to  first  approximation  as 
compared  to  that  from  X-rays  and  neutrons. 

For  the  NEO  materials  considered,  the  exponential 
model  indicates  that  when  D  AV  =  1, 
F(MT)Af2/Z)(km)2  is  smaller  for  X-rays  than  for 
neutrons,  while  for  D AV  =  10,  the  reverse  is  true  with 
neutrons  having  a  significant  advantage  over  X-rays 
in  this  range.  As  shown  in  Figs  2  and  3,  neutrons  are 
more  effective  than  X-rays  in  imparting  velocity 
changes  to  large  Fe/Ni  objects.  The  greater  pen¬ 
etration  of  the  neutrons  implies  that  more  mass  is 
ablated  and  hence  larger  yields  are  required  than  for 
X-ray  ablation,  however,  for  large  objects  (^  10  km), 
neutrons  would  appear  to  be  required  for  velocity 
increments  in  excess  of  a  few  cm/s.  The  exponential 
model,  for  a  given  diameter  and  D AV  product, 
requires  a  substantially  higher  yield  than  that  for 
constant  absorption  when  applied  to  ablation  by 
neutrons. 

From  eqns  (23)  and  (27),  the  maximum  impulse 
coupling  coefficient  occurs  for  FJ/16.6  =  0.386, 


y  =  1.166.  For  neutron  interactions,  this  implies 
DAV  «  280,  while  for  X-ray  ablation  of  VF  chon¬ 
drites,  NVF  chondrites  and  Fe/Ni,  one  has 
DAV  «0.4,  0.2,  and  0.04,  respectively. 

Due  to  the  difference  in  velocities  there  will  be  a 
tendency  for  the  X-rays  to  “fire  polish”  the  asteroid’s 
surface  before  arrival  of  the  neutrons.  Scattering  of 
the  neutrons  in  this  “atmosphere”  has  not  been 
considered  in  this  paper,  though  it  is  not  expected  to 
be  important  due  to  the  short  time  involved. 

The  above  conclusions  depend  critically  on  the 
assumed  radiation  spectrum  and  the  values  used  for 
the  opacities  as  well  as  the  very  simple  model.  Hence, 
the  results  should  be  taken  as  indicative  pending 
more  elaborate  estimates.  A  first-order  improvement, 
still  using  the  simple  model,  would  be  to  use  radiative 
opacities  averaged  over  the  Planck  distribution.  One 
would  then  substitute,  for  the  opacities  used  here, 
those  derived  by  suitably  averaging  them  using  the 
Rosseland  mean.  The  Rosseland  opacity  is  given 
by  [38]: 

Hr  =  (15/4jt4)| 

u  =  E/&  (31) 

where  E  is  the  X-ray  photon  energy  and  0  the  black 
body  temperature.  Thus  in  Section  3  one  would 
replace  n  by  fiR  wherever  it  appears.  If  we  further 
assume  that  n  ~  1  jEm,  then  it  is  easily  shown  that 
/iR~l/0m  and  hence,  using  eqn  (16),  the  specific 
impulse  should  vary  as  /  ~  0”"2.  For  a  hydrogenic 
plasma  m  =  3.5,  while  for  cold  matter  m  -  3. 
Integrating  over  all  the  absorption  lines  of  a  given 
material  can  lead  to  a  different  value  of  m.  For 
example,  from  Fig.  1 1  of  Shafer  et  al.  [20],  in  the  large 
fluence  limit,  m  =  2.0. 

Confirmation  from  a  subsidiary  calculation  as  to 
the  range  of  validity  of  the  exponential  model  in 
treating  the  neutron  effects  is  also  necessary  since  the 
sensitivity  of  the  neutron  results  in  the  simple  model 
used  herein  implies  that  more  elaborate  calculations 
are  desirable.  In  inertial  confinement  fusion,  the 
thermonuclear  energy  comes  from  a  highly  com¬ 
pressed  and  heated  target  consisting  of  a  DT  fuel 
surrounded  by  a  container.  Magelssen  and  Moses  [39] 
have  studied,  in  the  context  of  fusion  reactor  designs, 
the  X-ray  and  neutron  spectra  emitted  from  such 
targets  where  the  shell  was  a  high-Z  material  (Hg). 
The  spectra  they  find  are  consistent  with  those  we 
have  used.  Our  results  are  also  consistent  with  the 
work  of  Shafer  et  al.  [20],  both  for  X-rays  and  for 
neutrons.  Their  neutron  results  also  include  two 
points  derived  from  a  hydro-code  calculation.  In  a 
recent  paper,  a  team  at  the  Ecole  Polytechnique  [40] 
have  performed  experiments  on  the  impulse  imparted 
to  a  metallic  target  by  soft  X-rays.  Their  results 
confirm  the  general  trends  of  the  model  discussed 
above  as  applied  to  X-rays.  However,  more  exper¬ 
iments  along  these  lines  are  needed. 
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Finally,  recent  ground-based  radar  observations 
indicate  that  some  near-Earth  asteroids  are  far  from 
spherical  [41],  a  complication  that  will  have  to  be 
taken  into  account  in  future  simulations. 

The  potentially  important  case  of  nuclear  radiation 
interaction  with  comets  has  not  been  treated  explic¬ 
itly  in  this  paper,  but  may  be  approximated  by  using 
the  materials  properties  of  carbonaceous  chondrites, 
which  some  believe  are  similar  in  composition  to 
dormant  or  volatile  depleted  comets.  This  will  be 
discussed  elsewhere  using  a  newly  developed  model 
for  comets  [42], 
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